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Abstract. 

Hehum-like ions provide the most important X-ray spectral diagnostics in high 
temperature fusion and astrophysical plasmas. We previously presented computed 
collision strengths for O VII including relativistic fine structure, levels up to the n — A 
complex and radiation damping of autoionizing resonances. We have extended this 
work to other He-like ions (N, Ne, Mg, Al, Si, S, Ca). The calculations are carried 
out using the Breit-Pauli R-matrix (BPRM) method with a 31-level eigenfunction 
expansion. Collision strengths for the principal lines important in X-ray plasma 
diagnostics, w, x, y and z, corresponding to the 4 transitions to the ground level 
ls2 (i5o) < — ls2p{^P^),ls2p{^P^),ls2p{^P[),ls2s{^Si), are explicitly shown. We 
find the effect of radiation damping to be significant for the forbidden transitions in 
heavier He-like ions, which should affect the diagnostic line ratios. We extrapolated 
the collision strengths to their values at infinite energy using the Burgess- Tully 
extrapolation technique. This is required to calculate the Maxwellian average collision 
strengths at high temperature. We show that the coupling between dipole allowed 
and inter-combination transitions affects increasingly the effective collision strengths 
for the n^So — n/^Pi transition as the charge of the ion increases. This clearly affects 
the treatment of the extrapolation toward the infinite energy point of the collision 
strength. This work is carried out as part of the Iron Project-RmaX Network. 
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1. Introduction 

Helium-like ions provide the most important X-ray spectral diagnostics in high 
temperature fusion and astrophysical plasmas. The new generation of X-Ray satellites 
such as the Chandra X-Ray Observatory and the X-Ray Multi- Mirror Mission-Newton 
provide high resolution spectra of different types of astronomical objects (e.g. Kaastra 
et al. 2000, Porquet and Dubau 2000, Porquet et al. 2001). The high sensitivity of 
these observatories and the high quality of the spectra they produce requires highly 
accurate atomic data for a precise interpretation. The aim of the Iron Project-R-matrix 
calculations for X-ray spectroscopy (IP-RmaX) is to calculate extended sets of accurate 
collision strengths and rate coefficients for all ions of importance in X-Ray diagnostics. 
Among previous works, the electron impact excitation of Helium-like ions was considered 
by Sampson et al. (1983) and Zhang & Sampson (1987). They used the Coulomb-Born 
approximation with exchange, intermediate coupling, and some resonance effects to 
obtain collision strengths for Helium-like ions, with atomic number Z spanning a large 
range of values (4 < Z < 74). The present work aims at generating a more complete 
datasct of high reliability for He-like ions, including all important effects for highly- 
charged ions such as relativistic effects, radiation damping, and resonances in higher 
complexes up to n = 4. 

The method and computations are summarized in section 2. Results for the 
collision strengths and important issues are discussed in section 3, and the present 
results for the effective (MaxweUian averaged) coUision strengths are compared with 
previous calculations. The main conclusions are given in section 4. 

2. Method and Computations 

The collisional calculation in the present work has been carried out using the Breit-Pauli 
R-matrix (BPRM) method as used in the Iron Project (IP) and utilized in a number of 
previous publications. The aims and methods of the IP are presented in Hummer et al. 
(1993). We briefly summarize the main features of the method and present calculations. 

In the coupled channel or close coupling (CC) approximation the wave function 
expansion, ^{E), for a total spin and angular symmetry SLn or Jtt, of the (N-l-1) 
electron system is represented in terms of the target ion states as: 

^{E)^AY,XA + Y.<^,^3i (1) 

i 3 

where Xi is the target ion wave function in a specific state SiL^Tri or level JjTTj, and di is the 
wave function for the (N-l-1)*'* electron in a channel labeled SiLi{Ji)ni kfii{SLn) [Jvr]; kf 
is the incident kinetic energy. In the second sum the $j 's are correlation wave functions 
of the (N-l-1) electron-system that (a) compensate for the orthogonality conditions 
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between the continuum and the bound orbitals, and (b) represent additional short- 
range correlations that are often of crucial importance in scattering and radiative CC 
calculations for each symmetry. The <l>j's are also referred to as "bound channels", as 
opposed to the continuum or "free" channels in the first sum over the target states. In 
the relativistic BPRM calculations the set of SLir are recoupled in an intermediate (pair) 
coupling scheme to obtain (e + ion) states with total Jtt, followed by diagonalization 
of the (N-l-l)-electron Hamiltonian. Details of the diagonalization and the R-matrix 
method are given in many previous works (e.g. Berrington et al. 1995). 

The target expansion for the close coupling calculations consists of 31 fine-structure 
levels arising from the 19 LS terms with principal quantum number n < 4. The target 
eigenfunctions were developed using the SUPERSTRUCTURE program (Eissner et al. 
1974) in a version due to Nussbaumer and Storey (1978). The full expansion consists in 
10 configurations up to n=4: Is^, ls2s, ls2p, ls3s, ls3p, ls3d, IsAs, IsAp, IsM and ls4/ 
The scaling factors have been obtained by first optimizing the Ai^ parameter through the 
minimization of the energy of the ground term. Then all the other scaling parameters 
have been optimized on the 19 terms considered here. The scaling factors in the Thomas- 
Fermi potential employed in SUPERSTRUCTURE can be obtained upon request from 
the first author. The optimization of the orbitals has been performed as follow: 

In order to estimate the quality of the target wavefunction expansion, we compare 
the energy levels with those from the National Institute for Standards and Technology 
(NIST). The agreement with the NIST values is found to be very good, within 0.05% for 
all levels. A better criterion for the accuracy of the wavefunctions is the accuracy of the 
oscillator strengths for transitions in the target ion . The oscillator strengths agree well 
within 10% (however, some of the values are given by NIST with an estimated accuracy 
of 30%). Another accuracy criterion is the level of agreement between the oscillator 
strengths in the length and the velocity formulations, which we also find to be a few 
percent for all transitions. 

The 19 LS terms arc recoupled in the relativistic BPRM calculations into the 
corresponding 31 fine-structure levels up to the n = 4 complex using the routine 
RECUPD that performs intermediate coupling operations including the one-body Breit- 
Pauli operators (Hummer et. al. 1993). The reconstructed target eigenfunctions and the 
resulting target energies reproduce to 10^"^ Ryd the results from SUPERSTRUCTURE, 
verifying that the algebraic operations have been carried out self-consistently and 
without loss of accuracy. The collision strengths have been calculated for electron 
energies < < 4 x E{4}Pi) Ryd. This wide energy range ensures a good coverage 
of the region where resonances up to the n = 4 complex are important, as well as the 
higher energy region where no resonance has been included (all channels are open) but 
where the background collision strengths still make a significant contribution to the 
Maxwellian averaged rate coefficient for electron temperatures of interest. However, it 
might not be sufficient for the calculation of collision rates at high temperature. Given 
the slow convergence of some transitions and the need for high partial wave at high 
temperature, it is unpractical to carry the calculations at higher energy. However it is 
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possible to extrapolate the collision strengths to the value at infinite electron energy 
using the procedure from Burgess and Tully (1993). 

The inner region R-matrix basis set included 40 orbitals per angular momentum. 
Because of the importance of the near threshold resonances in the Maxwellian average 
rate coefficient, careful attention has been devoted to the resolution and a precise mesh 
has been chosen. A mesh of 10~^ Ryd was selected for the region where resonances 
are important, and a coarser mesh for the region where all channels are open. We 
included the contribution to the collision strengths from all symmetries with total 
angular momentum J and both odd and even parities, Jtt < (^)°'^. The contribution of 
higher partial waves was included using the Coulomb-Bethe approximation via the 'top- 
up' facility in the asymptotic region program STGF of the R-matrix package (Burke 
and Seaton 1986; modified by W. Eissner and G.X. Chen). 

3. Results and discussion 

In figures 1 and 2 we present the coUision strengths for transitions from 2^5*0 and 2^P" 
respectively to the ground state for all ions included in the present calculation. The 
high resolution of the calculations with a large number of points allows us to resolve 
clearly all the resonances up to the last threshold in the n = 4 complex. Wc delineate 
the Rydberg scries converging to the different scries limits in all three complexes. The 
doubly excited (c + He-like ion ) Li-likc ion resonance complexes, KMM, KMN etc.. 
converging towards the different n = 3 and n = 4 levels are clearly resolved. 

3.1. Radiation Damping 

It has been previously shown (Presnyakov and Urnov 1979, Pradhan 1981, Pradhan and 
Seaton 1985), that radiation damping may have a significant effect on the resonances 
in coUision strengths for highly charged ions since the radiative decay rates are large 
and compete with autoionization rates, i.e. the effect of dielectronic recombination 
on electron impact excitation. We studied in detail the radiation damping effect 
of dielectronic recombination, on resonance structures, collision strengths, and rate 
coefficients. We presented the detailed results for O vil in Delahaye & Pradhan (2002) 
(figures 3 to 6). We found for this ion that the radiation damping affected only a small 
region just below the threshold of convergence. As a consequence the effective collision 
strengths T are not affected at all. However, radiation damping is important for higher- 
Z elements since the transition probabilities increase with Z. Pradhan (1983a,b) has 
estimated the effect on the z-transition to be 9% in T for Fe xxv at the temperature of 
maximum abundance of helium-like iron, as also confirmed by the detailed calculation 
done by Whiteford et al. 2001. 

We calculated the effect of radiation damping and found it to be significant to more 
than 10% at low temperatures for several transitions between levels of Ca xix. As the 
temperature increases the effect on the rates decreases since the contribution of the 
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Figure 1. Collision strengths for the principal X-Ray line transition z (from ground 
state Is^ ^So to 2^Si). 



resonances near threshold to the rate becomes less important. 

In Figure 3 we present the details of the radiation damping effect on the collision 
strengths for the transition 2^5"! — 1^5'o of Caxix. We clearly see that the effect is 
important in the resonances near threshold and that it will affect the effective collision 
strengths at low temperature as we show in Figure 4. 
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Figure 2. Collision strengths for transitions corresponding to the w (or r) line 
(ls2 ^ 2ipf ) 



3.2. Effective Collision Strengths 

The Maxwellian averaged collision strengths 

T{T) = / n,,{e,)e-^^f''^d{e,/kT), (2) 
Jo 

have been computed for all transitions among levels up to the n = A. The collisional 
rates deduced for levels among n = 4 might not be as accurate as for level of lower 
complexes. Indeed, as has been shown for He-hke Fe XXV (Kimura et al. 1999, 2000, 
and Machado-Pelaez et. al. 2001), and for He-like O VII (Delahaye & Pradhan, 2002) 
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Figure 3. Collision strengths for the z-line (transition 2'^5i — 1^5o) in Ca XIX. The 
black line represents the present results when the radiation damping is neglected and 
the red one includes the effect. A zoom has been made on the region near threshold. 

the resonances arising from the complex n = N + 1 have a strong effect on transitions to 
and within the complex n = N. The extension to the n = 5 complex has been discussed 
by Whiteford et al. (2001) and will not be repeated. 

In Figure 4 we present the results of the damping effect on transitions 
(2^P°, 2^5"!, 2^P°, 2^P2, corresponding to the 4 principal lines w, z, x and y, that 
are of primary interest in X-ray spectral diagnostics (e.g. Gabriel and Jordan 1969, 
Pradhan 1982, Porquet et. al. 2001) . We can anticipate that the damping effect 
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will not affect significantly the density sensitive ratio R(Ne) = z/(x+y) because the 3 
transitions are affected similarly. However the ratio G(Te) = (x+y+z)/w, sensitive to 
both the temperature and the ionization balance in the plasma, will be reduced due to 
the unaffected w line compared to the 3 other affected lines. Therefore this temperature 
and ionization diagnostics would be significantly modified. 

We also compared the effective collision strengths with previous calculations from 
Zhang & Sampson (1987). The agreement is good with a tendency towards lower values 
in Zhang & Sampson and more pronounced for the z line. Indeed, the difference reaches 
a maximum of 17% at low temperatures for the z line, and 13% for the w line at high 
temperature. For the x and y lines the differences are below 10% and most of the time 
within 5%. 

From Figure 1 and 2 we can see that the general behavior is as expected. We can 
see Q decreasing as Z increases. It roughly follows 

Ooc(Z + l)-2. (3) 

This comes from the fact that in an isoclcctronic sequence the target is getting more 
compact as Z increases, the attraction from the highly charged nucleus being stronger 
for higher Z. We also sec how the complexes get farther apart as Z increases, decreasing 
the interaction between levels of the same symmetry but from different shells. 

In order to observe the general behavior along the isoelectronic sequence we use the 
reduced effective collision strengths as introduced by Burgess & TuUy (1992). It allows 
us to represent on the same plot the effective collision strengths for each transition for 
all ions considered in the sequence in a reduced temperature Tr — T/Z^ domain. 

The behavior at high energy normally should follow the behavior corresponding to 
the type of transition (Burgess & TuUy 1992). However, as we increase the charge the 
spin-orbit coupling is increasing and as previously shown (Pradhan 1983, also mentioned 
in Whiteford et al. 2001) in this case the rclativistic effects are crucial in order to account 
for this coupling, which of course docs not appear in the LS calculation. In the present 
case the coupling between the triplet P" and the singlet P" for example is becoming 
more and more important as we go to heavier elements. In the complex n — 2, the 
levels ^P° and ^P° are more mixed for heavier He-like ions. As a consequence the 
effective coUision strengths for transitions between the ground state and 2^P° which is 
an inter-combination transition should behave at high energy like 

n ~ const./ (4) 

and lim£;_,oo = but the mixing between the two levels, singlet and triplet, changes the 
behavior of this inter-combination transition by introducing a dipole-allowed component. 
The coulomb hmit then is not the same any more and does not tend to when E ^ 0. 
While Pradhan (1983) showed the trend for He-hke Fe, Se and Mo (Z=26,34 and 42 
respectively) the present work demonstrates that the rclativistic effects are important 
already at lower Z and not only for species with Z > 20 (Pradhan 1983). In Figure 5 
we can see the change in the behavior of the collision strengths as E increases for this 
inter-combination transition between the singlet S and triplet P°. While He-like N, O 
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Figure 4. Effective collision strengths for the principal lines (z, x, y, w). solid line: 
Present work with damping, Dashes: Present work without damping, dotted line: 
Zhang & Sampson 1987. 

and Ne present a discernible trend corresponding to inter-combination transitions, it is 
clearly obvious in He-like Si, S and Ca. This is of importance in the procedure to treat 
the infinite limit of the collision strengths. Indeed, the definition of the effective collision 
strengths (eq. 3.3) supposes an integration up to infinite energy. But in practice we 
have to stop at a certain point and make an extrapolation to infinity. As presented in 
Burgess and TuUy (1992), this extrapolation depends on the nature of the transition. 
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Figure 5. Detailed of reduced collision strengths for the inter-combination transition 
1^5*0 — 2^Pi. The triangles represent the Coulomb-Born limits as ^ oo. 

In Figure 6 we present a comparison of the effective collision strengths for the 
present results with 2 different treatments for the high energy limit for the transition 
1 ^5*0 — 2 ^Pi, as well as the results from Zhang & Sampson (1987) (in red). The solid 
lines assume a dipole behavior for the inter-comnination transition while the dotted lines 
indicate an inter-combination behavior without the relativistic mixing with the dipole 
allowed transtion l^^o - 2^P°. We can see that the former treatment is not adapted for 
the lighter elements and gives an overestimation of the effective collision strengths at 
high temperatures but is necessary for the heavier ones to avoid any underestimation 
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Figure 6. Effective collision strengths for He-like ions for transitions among the 
ground state and n = 2 complex. The solid black line assumes a dipole behavior 
at high energies for the inter-combination transition 1 — 2 '^Pf while the green 
dotted line is without relativistic mixing. In red are the results from Zhang & Sampson 
(1987). 



which would then directly affect the coUisional rates. 
4. Conclusion 

Some of the general conclusions of the paper are as follows. 
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1. The most complete close- coupling calculation using the Breit-Pauli R-Matrix 
method has been carried out for helium-like ions, including resonances up to n = 4 
levels. Detailed studies of radiation damping indicate that it may have a significant 
effect on the detailed collision strengths in a small energy region below the threshold(s) 
of convergence, but not on the effective collision strengths for the low Z ions of the 
sequence. However, radiation damping is important for higher-Z elements since the 
transition probabilities increase with Z and must be taken into account, especially at 
low temperature. 

2. The new results for the important X-Ray line transition may significantly 
affect the analysis of X-ray spectra from photoionized sources (e.g. active galactic 
nuclei), where ions may be abundant at relatively low temperatures. In coUisional 
ionized (coronal) sources the new results may not affect the theoretically computed line 
intensities significantly at temperatures close to maximum abundance, but should still 
do so at lower temperatures. It would be preferable to employ the present data in future 
coUisional-radiative and photoionization models. 

3. As all relevant atomic effects in electron-ion collisions have been considered, and 
resonances have been carefully delineated, we should expect the present results to be of 
definitive accuracy. Nonetheless, we conservatively estimate the precision to be about 
10-15%. 

4. A study of the line ratios R and G in He-like ions is in progress for diagnostics 
of astrophysical and laboratory sources. 

5. All data will be electronically available from the first author from 
franck.delahaye@obspm.fr. 

The authors would like to thank Dr. Werner Eissner and Dr. Peter Storey for 
their help and fuitful discussions. This work was supported partially by the NASA 
Astrophysical Theory Program. The computational work was carried out at the Ohio 
Supercomputer Center in Columbus, Ohio. 
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